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Leiden 

We study the sandpile model in infinite volume on 1,''. In partic- 
ular, we are interested in the question whether or not initial config- 
urations, chosen according to a stationary measure fi, are /x-almost 
surely stabilizable. We prove that stabilizability does not depend on 
the particular procedure of stabilization we adopt. In d = 1 and fj, a 
product measure with density p = 1 (the known critical value for sta- 
bilizability in ci = 1) with a positive density of empty sites, we prove 
that fj. is not stabilizable. 

Furthermore, we study, for values of p such that fi is stabilizable, 
percolation of toppled sites. We find that for p > small enough, 
there is a subcritical regime where the distribution of a cluster of 
toppled sites has an exponential tail, as is the case in the subcritical 
regime for ordinary percolation. 

1. Introduction. The sandpile model was originally introduced as a dy- 
namical model to illustrate the concept of self-organized criticality [1]. The 
model is defined on a finite subset A of Z*^, in discrete time. It starts with a 
stable configuration, that is, every site has a nonnegative height of at most 
2d—l sand grains. Every discrete time step, an addition of one sand grain is 
made to a random site. If this site becomes unstable, that is, has at least 2d 
grains, it topples, that is, it gives one grain to each neighbor. This may cause 
other sites to become unstable, and the topplings continue until every site is 
stable again. The total of all necessary topplings is called an avalanche, so 
after the avalanche we have reached the new configuration. This is possible 
in a finite number of topplings because, at the boundary of A, grains are 
dissipated. This model is abelian: the obtained configuration is independent 
of the order of topplings. 
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This sandpile model is said to exhibit self-organized critical behavior, for 
the following reasons. As the model evolves in time, it reaches a stationary 
state that is characterized, in the large-volume limit, by long-range height 
correlations and power law statistics for avalanche sizes, and thus reminds 
one of critical behavior in statistical mechanical models. However, the sand- 
pile model evolves naturally toward this critical state, without apparent 
tuning of any parameters. 

This seeming contrast has been discussed in [3, 4, 8]; it is argued that the 
model definition in fact does involve tuning. Namely, the instantaneousness 
of topplings, and the vanishing of dissipation as A | Z'^, can be viewed as 
a tuning of the addition and dissipation rate to respectively. This tun- 
ing then would ensure that the model evolves toward the critical point of 
a parametrized, nondynamical sandpile model, which can informally be de- 
scribed as follows. We start with an initial height configuration on Z*^ (not 
necessarily stable) according to a translation invariant probability measure 
with density p, which is the expected height, or number of sand grains 
per site. We keep toppling until there are no more unstable sites. If this 
is possible with a finite number of topplings per site, then we obtain the 
final configuration, and the initial configuration is said to be stabilizable. 
This version of the sandpile model was introduced in [3], and mathemati- 
cally investigated in [4, 8]. Results so far obtained are as follows: for d=l, 
any translation invariant probability measure with density p < 1 is stabiliz- 
able, any translation invariant probability measure with density p > 1 is not 
stabilizable and, for p = l, there are cases of stabilizability and nonstabiliz- 
ability; see [8]. For general d, any translation invariant probability measure 
with density p < d is stabilizable, and any translation invariant probability 
measure with density p> 2d — 1 is not stabilizable, in between d and 2d — 1 
there are nonstabilizable and stabilizable cases [4]. 

The present paper continues this investigation and from here on, when we 
talk about the sandpile model, we mean the version in infinite volume. In 
Section 2 we introduce notation, introduce general toppling procedures and 
discuss stabilizability issues. In this section we also prove that if a random 
initial configuration is stabilizable, then the expected height (density) is con- 
served by stabilization. In Section 3 we define critical values, and investigate 
the behavior at the critical point in d = l. We find that configurations cho- 
sen according to a nondegenerate product measure are a.s. not stabilizable. 
In Section 4 we investigate phase transitions for the sandpile model from 
a new viewpoint: we consider, for stabilizable configurations, percolation of 
the collection of toppled sites. We look for a critical p, not necessarily equal 
to pc mentioned above, such that, for all p below this value, there is no 
infinite cluster of toppled sites. 

For a general class of initial distributions and p small enough, we find 
a subcritical regime where not only there is a.s. no infinite cluster of top- 
pled sites, but the distribution of the cluster size has an exponential tail. 
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This corresponds to the subcritical regime for ordinary percolation, thus 
strengthening the idea of a critical phase transition. 

2. Toppling procedures and stabilizability. Denote hy X = N^^ the set 

of all height configurations and hj Q = {0, 1, . . . ,2d — 1}^ the set of stable 
height configurations. 

A toppling at site x applied to the configuration rj £ X is denoted by Ox (f?) 
and defined via 

f r]{y) - 2d, if y = x, 

(1) Ox{r]){y) = U{y) + l, if\y-x\ = l, 

[ r]{y), otherwise. 

A toppling at site x E Z"^ is called legal for the configuration t] £ X, if it is 
applied to an unstable site, that is, if r]{x) > 2d. 

The above definition of a toppling gives rise to the definition of the top- 
pling matrix A associated to the sandpile model. This is a matrix indexed 
by sites x,y £l/, with entries 

^x,y — 2d\x=y y|=l- 

With this definition, and with 5x defined to be the vector with entry 1 at x 
and entry in all other positions, we can write 

^xiv) =r]- A5x. 

Definition 2.1. A toppling procedure is a measurable map (with re- 
spect to the usual Borel sigma-algebra's) 

(2) T: [0,oo) X Z"^ X Af^N 
such that, for all i] £ X , 

(a) for all x G Z'^, 

T{0,x,ri) = 0. 

(b) for all X £ Z'^, 

t T{t, X, rf) 

is right-continuous and nondecreasing with jumps of size at most one, that 
is, for all t > 0, x £ U^, r] £ X, we have 

T{t,x,r]) - T{t-,x,r]) < 1. 

(c) for all X £'L'^ , in every finite time interval, there are finitely many 
jumps at X. 
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(d) T does not contain an "infinite backward chain of toppUngs," that is, 
there is no infinite chain of topphngs at sites Xj, i = 1,2, . . . , occurring at 
times ti > > • • • , where for all i, Xj+i is a neighbor of Xj. 

Note that condition (d) is only relevant in continuous time. We interpret 
T{t,x,rj) as the number of topplings at site x in the time interval [0,t], 
when T is applied to the initial configuration rj G X. The vector of all such 
numbers at time t is denoted by T{t,-,rj). We say that for all t such that 
T{t—,x,r]) < T{t,x,7]), site x topples at time t. 

If T is a toppling procedure, then for rj G X, t>0, we call 

e^(r) = {x e Z'^ : T{t, X, r]) > T{t-,x, r])} 
the set of sites that topple at time t>0 (for initial configuration rj). 

Definition 2.2. Let T be a toppling procedure. The configuration rjt 
at time t > associated to T and initial configuration r] £ X is defined to be 

(3) rit = i]-AT{t,;v)- 

Definition 2.3. A toppling procedure T is called legal if for all ij e X, 
for alH > and for ah x G ej'(r), ?7t-(x) > 2d. 

In words, this means that in a legal toppling procedure, only unstable sites 
are toppled. 

Definition 2.4. (a) A toppling procedure T is called finite for initial 
configuration rj G X, if for all x G Z'^, 

(4) T {oo , X , r]) := lim r(t, x, 77) = supT(t, x, r/) 
is finite. 

(b) A legal toppling procedure T is called stabilizing for initial configu- 
ration G ^ if it is finite and if the limit configuration 7700; defined by 

(5) r^oo =??- Ar(oo,-,r7), 
is stable. 

A random toppling procedure is a random variable with values in the set 
of toppling procedures. This can also be viewed as a measurable map 

T:[o,oo) xz'^ X X xn^n, 

where Cl denotes a probability space, and such that for all a; G except a 
set of measure 0, T(-, •, -,0;) is a toppling procedure. 



STABILIZABILITY IN SANDPILES 



5 



Definition 2.5. A toppling procedure is called stationary if for all t, the 
distribution of T{t, - ^rj) is translation invariant when we choose r/ according 
to a translation invariant probability measure. 

Next we discuss some examples. These examples have in common that 
for every t, if "qt contains unstable sites, then each of these sites will top- 
ple within finite time almost surely. As a consequence, for every 77, these 
toppling procedures are either stabilizing or infinite. Moreover, if they are 
infinite, then T{oo,x,rj) = 00 for every x. This can be seen as follows: if 
there is one site x that topples infinitely many times, then the neighbors 
of X receive infinitely many sand grains. Therefore, these neighbors need to 
topple infinitely many times, etc. 

1. Markov toppling processes. These are examples of random stationary top- 
pling procedures and are defined as follows. Each site rr S Z*^ has a Poisson 
clock (different clocks are independent) with rate one. When the clock 
at site X rings at time t and in the configuration rjt- , x is unstable, then 
X is toppled. More formally, the configuration r]t of (3) is evolving ac- 
cording to the Markov process with generator, defined on local functions 

L/(7?)=^l,(,)>2,(/(M)-/(r/)). 

X 

It is not hard to see that this procedure satisfies all requirements, in 
particular, (d) of Definition 2.1. 

It is also possible to adapt the rate at which unstable sites are toppled 
according to their height. In that case the Markov process becomes 

LJiv) = hi^)>2dc{v{x))if{0xri) - fiv)), 

X 

where c : N ^ M has to satisfy certain conditions in order to make the 
process well-defined. 

2. Toppling in nested volumes. This is a deterministic, discrete time toppling 
procedure. Choose a sequence Vn C Vn+i C Z'^ such that Un K = Z'^, but 
all Vn contain finitely many sites. We start toppling all the unstable sites 
in Vo until the configuration in Vq has no unstable sites left, then we do 
the same with Vi, etc. We put this into the framework of Definition 2.1 
as follows. At time t = 1, we topple all the unstable sites in Vq once, at 
time t = 2, we topple all the unstable sites in Vq if there are still unstable 
sites left after the topplings at time t = l, etc., until at time t = t{VQ,r]), 
no unstable sites are left in Vb; we then start toppling at time t = t{Vo,r]) 
all unstable sites in Vi, etc. Since the volumes Vn are finite, all t{Vn,ri) 
are finite. 
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We will use this procedure several times, but for ease of notation we 
will reparametrize time such that Vn is stabilized at time n instead of at 
time t(T4,?7). 

3. Topplings in parallel. Topplings in parallel consists simply in toppling at 
time t all unstable sites of r]t^i once. This toppling procedure is discrete 
time, deterministic and stationary. 

4. Topplings in waves. This procedure is only used for initial configurations 
having a single unstable site, say, at x^TL'^ . For the formal definition of 
the toppling function T{t,x,r]), we put it equal to zero for configurations 
rj that contain more than one unstable site. Toppling in waves is defined 
for the sandpile model on a finite grid as follows [6] : at t = l, we topple 
X once and, subsequently, all other sites that become unstable. All these 
topplings form the first wave. If after these topplings, x is still unstable, 
then at t = 2 we perform the second wave, etc. In each wave, no site 
topples more than once. 

This does not fit into our framework, because in each wave all topplings, 
except the toppling at x, are illegal. Nevertheless, we want each wave to 
be completed in finite time. Therefore, we define topplings in waves as 
follows: At t = 1, we topple site x once. Then for i = 2,3, . . . , at times 
2 — i, we consecutively topple all sites that are unstable except site x. 
That way the first wave is completed at time t = 2. All other waves 
proceed similarly. Since in each wave no site topples more than once, this 
procedure is well defined. 

Definition 2.6. (a) A configuration -q ^ X \s called stahilizahle if there 
exists a stabilizing legal toppling procedure. 

(b) A probability measure /i on [X,T) is called stabilizable if ^-almost 
every rj is stabilizable. 

Example 2.7. We give an example of a configuration that is not sta- 
bilizable: Consider the configuration ^ in Z where all sites have height 1, 
except the origin, which has height 2. From trying out by hand, it should 
become clear that this configuration is not stabilizable. We may choose to 
topple in waves, since there is only one unstable site. In our case, in each 
wave every site topples exactly once, so that after each wave we obtain the 
same configuration; hence, there are infinitely many waves. Alternatively, we 
may choose to topple in parallel. Then in our case, the height of the origin 
alternates between and 2, so that the origin topples infinitely many times. 
From the forthcoming Theorem 2.8, we can use either (3) or (4) to conclude 
that ^ is not stabilizable. 

In [4], Definition 2.4, stabilizability is defined in terms of toppling in 
nested volumes, and in [4], Lemma 6.12, it is proved that this definition 
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of stabilizability is equivalent for this toppling procedure and the Markov 
toppling procedure. Here, we extend this result: we prove that if rj is stabi- 
lizable, then irrespective of what legal toppling procedure we choose, it will 
always be finite, and irrespective of what stabilizing procedure we choose, 
we always obtain the same stable configuration. On the other hand, if we 
find one infinite legal toppling procedure for r], then we know that 7] is not 
stabilizable. This can also be concluded from the existence of a legal toppling 
procedure in which every site topples at least once. 

Let T, T' be two toppling procedures, which are finite for initial configu- 
ration rj. Then we write 

T' ^r, T, 

if for all X G Z'^ 

T'(oo, X, rj) < T{oo, X, rf). 



Theorem 2.8. Let T , T' he two legal toppling procedures, which are both 
finite for initial configuration rj: 

1. IfT is stabilizing for Tj, then 

(6) r T. 

2. If T and T' are two stabilizing toppling procedures for rj, then for all 
x^T^, 

r'(oo, X, T]) = T(oo, X, rf). 

In particular, this means that for stabilizable r], the limit configuration 
rjoo is well defined. 

3. For stabilizable r] £ X , there does not exist a nonfinite legal toppling pro- 
cedure. 

4. // T is stabilizing for r], then there is at least one site x that does not 
topple, that is, there is at least one site x for which T[oo,x,rj) = 0. 

Proof. The proof of Statement 1 is inspired by an argument that ap- 
pears in [2] and [9] in the context of finite grids or discrete time toppling 
procedures. 

For every x, we define a time := sup{t : r'(t, x, 77) < T(cxd, x, r;)}, and 
we call all topplings in T' that occur at times strictly larger than Tx "extra" 
topplings. We suppose the converse of Statement 1, that is, we suppose that 
there is at least one extra toppling. 

Suppose an extra toppling occurs at site y, at time ty < 00. Then just 
before time ty, the number of topplings at site y is at least T{oo,y,rj). 
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Moreover, in order for this extra toppling to be legal, site y must be unstable 
just before time ty. Thus, we find, following T', that 

2d < r^ty^y) = 7]{y) - {AT'{ty-, •, r?))(y) 

= viy) - 2dT'{ty-,y,r]) + '^T'{ty-,x,r]) 

< r]{y) - 2dT{oo,y,r]) + ^^T' {ty-,x,r]), 

x^y 

where the sum runs over all neighbors of y. Since T is stabilizing, we 

have 

2d > r]{y) - 2dT{oo, y,v) + Yl ^' 

xr^y 

so for at least one x ~ y, T'{ty—,x,T]) > T{oo,x,t]). In other words, for an 
extra toppling at site s to be legal, it is necessary that it is preceded by at 
least one extra toppling at one of its neighbors. Then for this extra toppling, 
we can make the same observation. Continuing this reasoning, we find that 
in order for the extra toppling at s to be legal, we need an infinite backward 
chain of extra topplings, occurring in finite time. But then T' docs not satisfy 
item (d) of Definition 2.1. This proves Statement 1. 

To prove Statement 2, we simply observe that if T and T' are both stabi- 
lizing, then according to the above, T' T and T ^^^T' , so that they must 
be equal. 

To prove Statement 3, let T be a stabilizing toppling procedure, and T" a 
nonfinite legal toppling procedure. Since T" is nonfinite, there exists x S Z'^ 
such that T"(t, x, r/) | cxo as 1 1 c«. For some w < oo, we define as follows: 
for all t<w, Tll,{t,-,ri) = T"{t,-,ri), but for alH > u;, Tll,{t,-,r]) =T"{w,-,ri). 
In words, performs all topplings according to T" up to time w, but 
then stops toppling. is a finite legal toppling procedure by item (c) of 
Definition 2.1 and, hence, by Statement 1 of this theorem, T!^{oo,x^rf) < 
T(oo,x,r/). By letting w — > oo, we obtain a contradiction. 

To prove Statement 4, suppose that there is a stabilizing toppling proce- 
dure T such that T(oo,x,r/) > for all x. For every x, we call the toppling 
that occurs according to T at time tx :=min{t : T{t,x,rj) = T{oo,x,r])} the 
"last" toppling. Since T is stabilizing, is finite for all x. 

We define T as 

T{t, X, 7]) := min{T(i, x, r]),T{oo, x, rj) — 1}, 

so that for all x, T{oo,x,rj) = T(oo,x,r]) — 1. In words, T contains all top- 
plings according to T except the last one at each site. Note that T is a finite, 
but not a priori legal toppling procedure. However, we have 

T] - Af{oo,-,r]) = r]- AT{oo,-,r]) =1]^, 
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so that after all topplings according to T, we have a stable configuration. 
Now the argument proceeds as in the proof of Statement 1: we have, for 
some site v, 

2d<i]t^_{v) =T]{v) -2dT{t^-,v,r]) + ^T{ty-,x,r]) 
= rj{v) -2dT{oo,v,'n) + ^ T(t„-, x, ry), 

whereas 

2d > r}{v) — 2dT{oo, v, rf) + ^ ^"(00, x, rj). 

x^v 

Similarly, as in the proof of Statement 1, we conclude that for the last top- 
pling at V to occur legally, it must have been preceded by an infinite back- 
ward chain of last topplings, so that T cannot satisfy item (d) of Definition 
2.1. □ 

Remark 2.9. Note that if ^ is stabilizable and ergodic, then the induced 
measure on limit configurations is also ergodic since it is a factor of fi. 

We now prove that a finite legal toppling procedure conserves the density. 
From here on, we will denote by E^, expectation resp. probability with 
respect to /x. 

Lemma 2.10. Let fj, be a translation invariant probability measure on 
X such that £^(77(0)) = p < 00. Suppose, furthermore, that p is stabilizable. 
Then the expected height is conserved by stabilization, that is, 

E^(r?oo(0))=p. 

Moreover, if p is a translation invariant probability measure on X such that 
E^(?/(0)) = cxD, then fi is not stabilizable. 

Proof. Using ergodic decomposition, we can assume that fi is ergodic. 
We start with the case E^(r/(0)) = p < 00. Without loss of generality, we 
assume that the toppling procedure that stabilizes n is stationary, and, is 
moreover, such that for all t, E,^(T{t,x,rj)) < 00 [we can, e.g., choose the 
Markov toppling procedure, where T{t,x,rj) is dominated by a Poisson pro- 
cess] . 

At time t, we then have 

rit{x) = r]{x) - l^x,yT{t, y, rj), 
y 

which upon integrating over the distribution of r/ gives 
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Therefore, using Fatou's lemma. 

Poo :=E^(7?oo(0)) = Ej lim ??t(0)') < liminf E^(r?t(0)) = p. 

The inequahty E^(?7oo(0)) > p is proved in [4]; we give a somewhat different 
argument here. Let Xn denote the position of a simple random walk starting 
at the origin (independent of r]), and denote E,rw,Prw expectation and prob- 
ability with respect to this random walk. We start by choosing a stabilizable 
r] with limit rjoo, and for a moment we consider this r] and ?/oo fixed. 
From the relation 

Voo{x) = 7]{x) - ^ A^,j^T(oo,y,?7), 
y 

we obtain 

1 /""^ \ 2d 2d 

-E™ V(r/oo(^fc)-??(Xfc)) =—Ern^{T{oo,Xn,r])) r(oo,0,7?). 

n J n n 

(7) 

By letting n — > oo, using that r(oo,0, ry) < cxd by assumption, and that 
Krw{T {oo , Xn, r])) > 0, this leads to 

(8) liminf -E,^ r]oo{Xk)] > liminf -E,^ (j2 vi^k)] ■ 

If we now finally choose t] according to p, which is ergodic, then the limiting 
measure is also ergodic according to Remark 2.9. By ergodicity of the scenery 
process {ri{Xn) :n G N} (see, e.g., [4], Proposition 8.1), it follows that for 
a.e. r], the right-hand side is equal to p, and the left-hand side is equal to 
Poo- This proves that poo ^ P- 

Finally, if E^(7y(0)) = oo, then the right-hand side of (8) goes to -|-oo as 
n ^ oo. Therefore, if p would be stabilizable, then rjoo would have infinite 
density, a contradiction. □ 

3. Criticality and critical behavior. Let V{X) denote the set of all trans- 
lation invariant probability measures on We say that a subset M of 
V{X) is density complete if for all p S [0,oo) there exists p £ M such that 
p(r?(0)) = p. 

Let Ai C V{X) be density complete. We define the M-critical density for 
stabilizability to be 

(9) Pc{M) =sup{p> 0:VpG with p{r]{0)) =p,p is stabilizable}. 

Of course, it can be questioned whether the density is the only relevant 
parameter distinguishing between stabilizability and nonstabilizability. It is 
certainly the most natural parameter, and is considered in the numerical 
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experiments of [3]. In [4] a related notion of maximal stabilizability is intro- 
duced. 

It is clear that pc{M) < Pc{M') for M 3 M' . Natural choices for M 
are a one-parameter family of product measures such as the set of Poisson 
product measures with parameter p, the set of all product measures or simply 
M=V{X). 

The following results are reformulations of results in [4] and [8]. 

Theorem 3.1. (a) For M = V^X), and for all d, 

Pc{M) = d. 
(b) For all M density complete, we have 

d<pc{M) <2d-l. 
In particular, when d = l and for all M. density complete, we have 

Pc{M) = 1. 

We now specialize to the case d = l. Accordingly, let be a one-dimensional 
translation invariant product measure with density p = 1- From Theorem 
3.1, we know that for p < 1, ^ is stabilizable, and that for p > 1, it is not. 
The next result deals with the critical case p= 1- 

Theorem 3.2. Let p be a one- dimensional product measure with p = l 
such that p{r]{0) = 0) > 0. Then p is not stabilizable. 

Our strategy will be to show that there a.s. exists a nonfinite legal toppling 
procedure. This implies, using Theorem 2.8, that p is not stabilizable. In 
order to do so, we will use topplings in nested volumes, but for the proof 
it will be important to define an intermediate toppling procedure, during 
which we only stabilize in volumes of the form [0,n], that is, we increase 
the stabilized volume only to one side. After stabilization of the interval 
[0,n], the outer boundary sites —1 and n + 1 possibly contain, on top of 
their original height, extra grains that were removed from the interval [0,n] 
during stabilization; all other sites outside [— l,n-|-l] still have their original 
height. 

For a while we concentrate on this one-sided procedure. In this section we 
denote by rjn the configuration that results from stabilizing in the interval 
[0,n]. As in the nested volumes toppling procedure, we re-define time as to 
match this notation: stabilization of [0,n] takes place at time n so that rjn 
is the configuration reached at time n. 

We will work with the number and positions of empty sites of rjn in [0, n] , 
and we will call such an empty site "a 0" of rjn ■ In Figure 1 we illustrate the 
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Fig. 1. The first 10,000 time steps for stabilizing 77 according to Poisson(l) product 
measure in nested volumes [0,?i]. The y-axis represents time. A black dot indicates an 
empty site (a "0"). In addition, the outer boundary sites of [0,n] are also colored black. 
See the text for further explanation. 

dynamics of the O's in this procedure. Time (in the new sense) is plotted 
vertically going downward; space is plotted horizontally. At every time, when 
you look horizontally, the black dots you see represent the positions of the 
O's at that time in the interval [0,n]. In addition, the outer boundary sites of 
[0,n] are also colored black. The configuration outside the stabilized interval 
is not shown. Thus, the picture does not give complete information about the 
configuration rjn, it only shows the positions of the O's and the width of the 
stabilized interval. One can clearly follow the trajectory of the leftmost 0. For 
instance, just before time 4000, the leftmost starts to move, and reaches 
the outer boundary somewhere between times 6000 and 7000. Shortly after 
that, a new leftmost starts to move to the right, etc. 

Our strategy is to show that during the one-sided procedure, despite the 
fact that infinitely often new O's are created, we infinitely often encounter a 
configuration that does not contain a 0. Every time this occurs, there is a 
fixed positive lower bound for the probability that the origin topples. This 
will then imply that the origin topples infinitely many times a.s. 

In order to show that infinitely often there are no O's, we need to analyze 
the dynamics of the O's in great detail. We are going to view the O's as 
objects that can move, disappear or be created. In order to precisely define 
these terms, we organize the topplings that occur in time step n into waves. 
If at time n — 1 site n is unstable (at time n — 1, all sites 0, . . . ,n — 1 are 
stable), then in the first wave at time n, we topple site n once and then 
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all other sites in [0,n] that become unstable, except site n again. If after 
this wave site n is still unstable, then the second wave starts, etc. We will 
number the waves A; = 0, . . . and call fjn-i^k the configuration after the 
kth wave, so that fjn-ifl = rjn-i and fjn-i^K = Vn- Depending on the position 
of the rightmost after wave A; — 1, wave k has the following effect: 

1. If the rightmost of ??n-i,fc-i is at site n — 1, and after wave k site n is 
not empty, then the number of O's has decreased by 1; 

2. If there are no O's in fin-i^k-i, then all sites in [0,n] topple, after which 
there is a at the origin, and site —1 has gained a grain. 

3. In all other cases, if the rightmost is at position x — 1 (so that x is the 
leftmost site that topples during wave k), then site x — 1 gains a grain 
and site x loses one. In addition, site n loses one grain and site n + 1 
gains one. 

These observations inspire the following definition. 

Definition 3.3. Let, at time n, K be the number of waves. If K > 0, 
then let, in wave k, x he the leftmost site that topples: 

• If X = n, and after wave k site n is not empty, we say that the at site 
n — 1 disappears. 

• If X = 0, we say that a new is created at the origin. 

• If X > 0, and no disappears, we say that the at site x — 1 moves to 
site X. 

If r]n-i{n) = (this implies K = 0), we say that a new is created at the 
right boundary. 

Since there may be multiple waves in one time step, multiple things can 
happen to the O's. However, note that we have the following restrictions: in 
each wave, only the rightmost can move. For instance, in the example in 
Figure 1, the that is present at position 472, at time 10,000, has been in 
that position for almost 2000 time steps, and we cannot be sure if it will ever 
move again some future time (actually, as our proof will show, it will a.s.). 
Furthermore, only when, after a previous wave, there are no O's left can a 
new be created at the origin. For instance, in the realization in Figure 1, 
this occurs seven times between n = 6000 and n = 8000. 

We stress that, according to the above definition, we actually identify 
certain O's in different time steps. A look at Figure 1 should convince the 
reader that this is a natural way to view the O's, even though in order to do 
so, it is necessary to break the topplings in each time step up into waves to 
ensure a correct identification. Once a has been created, it exists until it 
disappears at the right boundary. This may be in the same time step, but 
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it could also require many time steps. During this time, it may move to the 
right or remain for some periods of time in the same position. 

The time intervals between successive instances where the number of O's 
is equal to some given number z are not i.i.d. time intervals. However, we 
will show in the following lemma that, for all z > 0, the time intervals, from 
the moment that, the number of O's becomes 2 + 1 until the first return to 
a value that is at most z, are i.i.d. time intervals, whose distribution does 
not depend on z. In the proof we use that, for z > 0, the number of O's can 
only increase from z to z + 1 when a new is created at the right boundary. 
When z = we can have that the number of O's increases because a new 
is created at the origin, in which case the proof does not apply. 

Lemma 3.4. Let z > 0. Let Z[n) be the number of O^s after time n. For 
i = 0, 1, . . . , let No{z) = 0, Mi{z) = mm{n > Ni{z) : Z{n) < z}, and Ni{z) = 
min{n > Mj_i (z) : Z(n) = z + 1} . Then: 

1. The random variables Aj(z) = Mi[z) — Ni{z) are i.i.d., for all i > 0. 

2. The distribution of Ai{z) does not depend on z, and we denote by A a 
random variable with this distribution. 

3. If liminf oo Z{n) < oo a.s., then liminin^oo Z{n) < 1 a.s. anc/ P(A = 
oo) = 0. 

Proof. Since at all times Ni{z) a new is created at the right boundary, 
it must be the case that 

miz)-i{Ni{z)) = o. 

This new will be present until time Mj(z) and, during this time, it cannot 
move to the left. The key observation is that the dynamics of this new 
depend only on r]{j), j > Ni[z). In particular, conditioned on the creation of 
the new at time Ni{z), the value of Mi{z) only depends on these random 
variables. Since 

1. the Mi{z) and Ni{z) are all stopping times, and 

2. is a product measure, 

it follows that the Aj(z) are i.i.d. random variables for every fixed z. Fur- 
thermore, it follows immediately that the distribution of Aj(z) is also inde- 
pendent of z. This proves the first two parts of the lemma. 

We now proceed with Part 3. Suppose that liminf^^oo Z{n) > 1 with pos- 
itive probability, that is, there is a random variable A^, finite with positive 
probability, such that Z{n) > 1 for all n> N. We denote by Ne{l) the total 
number of time intervals Aj(l) [during which Z{n) > 1]. If Ne{l) < oo, then 

the last interval has infinite length. However, sinc6 IP(lim inf^ Z{n) > 1) > 

0, it is the case that P(A'^e(l) < oo) > 0. We calculate, using that the Aj(l)'s 
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are independent, F{N,{1) < oo) = ^^=1 HNeil) = k) = ^^=1 nf=i' HM^) < 
cx))P(Afc(l) = oo). This equals J2'^^-^F{A < oo)''-1F(A = oo), so that we ob- 
tain P(A = oo) > 0. 

So far, we showed that hminf„__>oo ■^(^) > li with positive probabiUty, im- 
phes that P(A = oo) > 0. Now we show that P(A = oo) > imphes 
Hminf„^oo Z{n) = oo a.s. Denote by Nf>{z) the total number of time inter- 
vals Ai{z) [during which Z{n) > z], and call P(A = oo) =p. Similarly to the 
above computation, we calculate ¥{Ne{z) < oo) = X^fclilPlA < oo)''-^P(A = 
oo) = Efclo(l -^5)^ = 1- Therefore, P(iVe(z) < oo) = 1 for all z > 0, so that 
liminf„_»oo Z{n) > z a.s. for all z > 0. It follows that liminf„_>oo Z{n) = oo 
a.s. □ 

Proof of Theorem 3.2. We choose ry according to fi. As mentioned 
before, we will show that there a.s. exists a nonfinite legal toppling proce- 
dure. We will use toppling in nested volumes [—m,m]. However, in order to 
compare this procedure with the one-sided procedure introduced above, we 
will reach from r] in the following way: First, we stabilize the interval 
[0,m]. After this step, site —1 received a number A'^{m) of grains, and [0,m] 
contains a number Z~^{m) of O's. Then, we stabilize the interval [— m, — 2], 
in the same way as we stabilized in [0,m]. After this step, site —1 received 
another number A~{m) of grains, and [— m, — 2] contains a number Z~{m) 
of O's. Site —1 is now the only possibly unstable site in [— m,m]. Finally, we 
stabilize all of [—m,m]. Note that in this description, for every m, we obtain 
r/m starting from i], whereas in the above presented one-sided procedure we 
obtained r?^ from The numbers A~^{m) and A~{m) are nondecreasing 

in m. The sequences (Z"*"(m)) and {Z~{m)) are independent of each other, 
and also have the same distribution. 

The following discussion will repeatedly involve both the one-sided and 
the nested volume toppling procedure. To make the distinction clear, we will 
use indices n or to refer to time steps for the one-sided procedure, and 
indices m or M to refer to time steps of the nested volume procedure. 

First case. We assume that, with positive probability, liminf„^oo Z^{m) = 
00. If both liminfs are actually infinite and we apply the right one-sided pro- 
cedure to 7], then for every z > there is a time N{z) such that, for all 
n > N{z), [0, n] contains at least z O's. This, however, implies that the left- 
most z — 1 O's never move again, which in turn implies that from some n 
on, grains can never reach site —1 again; a similar argument is valid for the 
left one-sided procedure on the interval [— n, — 2]. Hence, there is positive 
probability that both A^{m) and A~{m) do not increase anymore eventu- 
ally, and therefore remain bounded. However, since both Z^[m) and Z~[m) 
tend to infinity, the number of O's with fixed positions in both the left and 
right one-sided procedure tends to infinity. This now is incompatible with 
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stabilization, since after toppling site —1 in the end, we should (if stabiliza- 
tion occurs) obtain a stable configuration r/oo which should be equal to I, 
by Lemma 2.10. However, there are simply not enough grains at —1 to fill 
all the O's that were created by the one-sided procedures. 

Second case. We now know that liminfm^oo Z^{m) < co a.s. By Lemma 
3.4, part 3, we conclude that a.s. liminf^^oo Z{n) < 1, and P(A = oo) = 0. 
This implies that all O's, possibly except the leftmost one, will eventually 
disappear. Although the proof of Lemma 3.4, Part 3, does not work when 
z = 0, we can a fortiori conclude that also the leftmost must eventually 
disappear. Indeed, since all other O's eventually disappear, and since the 
occurrence of this event only depends on the configurations to the right of 
such a 0, it follows that no matter what the configuration to the right of a 
certain is, it will always disappear eventually. Clearly, this is then also true 
for the leftmost 0. [Note though that the leftmost may disappear without 
Z{n) decreasing; if in a time step where the leftmost disappears also the 
origin topples, then a new is created at the origin.] 

Finally, since clearly infinitely many O's are created at the right boundary, 
we conclude that infinitely often the leftmost disappears. Now consider one 
time instant A'"' such that the leftmost disappears at time N' . Whether 
a new is created at this time depends on the precise value of 7?7V'-i(-^')- 
Given that this amount is large enough to make the leftmost disappear, 
we can either have that the origin topples as well, or we can have that the 
origin does not topple, so that at time N' there are no O's. In the last case, 
if r]]\f'{N' + 1) > 2, then the origin topples at time A^' -|- 1. The probability 
that 7]N>{N' -M) > 2 is bounded from below by P(7?(A' + 1)>2). Thus, we 
have that at every time instant where the leftmost disappears, either the 
origin topples, or it topples with at least a fixed positive probability one 
time step later. We conclude that during the one-sided procedure the origin 
topples infinitely often, so that the procedure is nonfinite. □ 

If the initial measure fi satisfies a central limit theorem, then a shorter 
proof is possible: 

Theorem 3.5. Let ^ he a translation invariant probability measure on 
X such that E^(r/(0)) = 1 and such that -^J2x=-ni'n{^) ~ 1) converges in 
distribution, as n ^ oo to a nondegenerate normal random variable. Then 
^ is not stabilizable. 

Proof. Suppose that fx is stabilizable. Then there exist random vari- 
ables N{x) such that, for all x G Z, 



(10) 



77 - AN{x) = 1, 
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where AN{x) ■.= J2y^xy^{y), and where 1 denotes the configuration with 
height 1 everywhere. Indeed, if is stabihzable, then the only final stable 
configuration can be the configuration which is constant and equal to 1 [by 
conservation of density (Lemma 2.10) and stability]. The variables N{x) are 
then the number of topplings needed at x to stabilize rj. By stationarity of 
the toppling mechanism, the joint distribution of N{x) is stationary under 
translations. From (10), we obtain 

Y: (r?(x) - 1) = ( Y: ^Nia 

(11) 



x=—n V \x=—n 



-^(Ni-n - 1) - Nin) + N(n + 1) - N(-n)). 



We claim that the right-hand side of (11) converges to in probability as 
n — > oo. Indeed, for e > 0, by stationarity, we have 



(-^(N(-n - 1) - N(n) + N(n + 1) - N(-n)) > e) 



<4P^(^iV(0)> 

which tends to zero by the assumption that A^(0) < oo a.s. This leads to a 
contradiction since, by assumption, the left-hand side of (11) converges to a 
nondegenerate normal random variable. □ 

4. Sandpile percolation. We call % the set of all sites that have toppled 
at least once up to (and including) time t, that is, Tt = {x: Tt{x) > 0}. Like- 
wise, we introduce the set of nonempty sites at time t, Vt = {x: r]t{x) > 0}, 
and finally Wj = TtDVt, the set of sites that have toppled or are nonempty 
at time t. 

For rj stabilizable, these sets have a limit, for example, 7^ = lim(__+oo 
We decompose the set 7^ in clusters '7^(x), where ^^(x) is the largest 
connected component of T^o containing x. Likewise, we decompose Woo into 
clusters Woo (a;), where Woo(ic) is the largest connected component of Woo 
containing x. Sandpile percolation is the study of these clusters. 

As in classical percolation, one can define critical densities for the ex- 
istence or absence of infinite clusters and distinguish between a sub- and 
supercritical regime. In this section we are interested in the tail of the clus- 
ter size distribution 

P,.(|Too(0)| >n) 
and in the percolation probability 

Pm(I^oo(0)| =oo). 
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For the other sets, definitions and notation are similar. In this section we 
will need the following large deviation result. 

Lemma 4.1. Let fip, < p <1 be a one-parameter family of translation 
invariant product measures satisfying the following: 



(i) For all t € M, the moment generating function 



(12) G,(t)=E^,(e*''W)<oo 



exists; 

(ii) The first moment satisfies 

IEm,(??(0))=p; 

(iii) If ap>0 are nonnegative such that Op ^ oo for p — > 0, then 

Gp^ (op) := sup{x : Gp{x) < Op} 

also tends to oo as 0. 

Then we have, for any sequence xi,X2, ■ ■ ■ of lattice points and any e > 0, 

1 / " \ 

(13) lim sup lim sup — log ^?/(xj) >en\= —oo. 

p^O n.-+oo n \i=l J 

Proof. By the Markov inequality, for any t > 0, using that under the 
r]{xi) are i.i.d., we obtain 



which gives 



J2v{x^)>en)<e-''''Gpitr, 



(14) ^logFp^\^^i]{xi)>enj < -et + logGp{t). 



n 



We now show that there exists tp> Q such that 

\ogGp{tp)<l 

and such that tp — > c« as p — > 0. Using that 7/(0) takes only integer values, 
the elementary inequality log(l + x) < x, the Cauchy-Schwarz inequality 
and, finally, the Markov inequality, we obtain 

logGp(t)<log(l+Ep^(e*''Wl,„>i)) 
<Ep^(e*''(°)l,(o)>i) 

(15) 

<Gp{2tfl^p^'\ 
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SO we can choose tp = ^Gp^{l/ p^^"^) and, by condition (iii), we have that 
tp — > oo as p — > 0. We can now finish the proof by returning to (14), choosing 
t = tp in the right-hand side of the inequaUty, and letting /) ^ 0. □ 

We remark that an elementary computation shows that conditions (i)-(iii) 
are satisfied when p,p has Poisson-p one-dimensional marginals. 
We now first deal with the tail of |7^(0)|. 

Theorem 4.2. (a) Let d=l and let p, he a translation invariant product 
measure satisfying Ep(e*''('^^) < oo and with density p <1. Then there exists 
a constant ci > such that 

PM(l^oo(0)|>?i)<e-=i". 

(b) Let d> 1 and let pp, < p < 1 be a collection of product measures 
satisfying conditions (i)-(iii) of Lemma 4-1- Then for all p sufficiently small, 
there exists a constant c^ = Cd{d,p) > such that 

PM,(|Too(0)|>n)<e-^''". 

For the proof, we need the following result which goes back to at least [7] . 
For completeness, we give the short proof. 

Lemma 4.3. Let A be a subset of Tt, for some toppling procedure. Let 
P\ be the number of internal bonds in A, that is, bonds with both endpoints 
in A. Then 

Proof. For each internal bond of A, consider the last particle that 
traversed this bond via a toppling. This particle remains in A up to time t. 
The result now follows. □ 

Proof of Theorem 4.2. Let /x be a product measure in dimension 
d>l. We then have 

oo 

P/.(|Too(0)| >n)= ^ Pp(|Too(0)| =m) + P^(|Too(0)| =oo). 

m=n 

We choose to stabilize in nested boxes -B^ of radius k. Recall that we 
reparametrize time so that, at time k, the whole box has been stabi- 
lized. Then for every k, the maximum size of Tk{0) is {2k + 1)'^, so that we 
can rewrite 

{2k+l)'' 

(16) Fp(|roo(0)| > n) = lim Pp(|rfc(0)| > n) = lim ^ PM(l'7'fe(0)| = m). 
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We will derive a bound for P^(|7^,.(0)| = m). We write 

\C\=m 
06C 

where the sum runs over all finite connected subsets C of size m containing 
the origin. Then, by Lemma 4.3, this implies a minimum number of at least 
m — 1 sand grains in C in r/oo • But since no sand can have entered C during 
stabilization — in fact, grains must have left C — it also implies that C contains 
at least m grains at t = 0. Since ^ is a product measure, this corresponds 
for p < 1 to a classical large deviation of ^Z)xeC^(^)' bound 
the corresponding probability by a Chernov bound for sums of independent 
random variables, that is, there is a constant a = a{p) > such that, for all 
m, 

(17) P^(rfc(0)=C)<e-"™. 

(For this statement we do not yet need Lemma 4.1.) We now distinguish 
between (a) and (b): 

(a) For d= 1, the number of clusters of size m containing the origin is 
equal to m. Hence, for d= 1, we arrive at 

{2k+l)'^ {2k+l)'^ {2k+l)'^ 

J2 F^(|rfc(0)|=m)< E^-""^^ E me-"- < e-'^^-, 

m=n m=n \C\=m m=n 

oec 

with ci positive for all p < 1. Since this outcome does not depend on k, when 
inserting this in (16), we obtain for d = l 



lPM(|Too(0)|>n)<e- 



proving (a). 

(b) For d > 1, we have, according to Lemma 4.1, that limp^Q a{p) = oo. 
Also, for d> 1, there is a constant a' = a'{d) such that the number of clusters 
of size m containing the origin is at most e" see, for example, [5]. Hence, 
we calculate 



{2k+l)'^ {2k+iy (2fc+l) 

E 



m=n m=n \C\=m m=n 

OGC 



with Cd positive for p small enough. The proof is now finished as in case (a). 
□ 
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Remark 4.4. In the proof of Theorem 4.2, we used that is translation 
invariant, and that we have, for p smah enough, a large deviation bound for 
sums like J^xeAVi^)^ with A some connected volume in Z*^. There are many 
more measures that satisfy these requirements, for instance, Gibbs measures 
or other sufficiently rapidly mixing measures. 

The argument to prove the exponential tail of the distribution of |yVoo(0)|, 
which in turn implies the exponential tail of the distribution of |Voo(0)|, is 
similar, although some extra arguments are needed. 

Theorem 4.5. Let fip, < p < 1 be a collection of product measures 
satisfying conditions (i)-(iii) of Lemma 4-1- Then for all p sufficiently small, 
there exists a constant = ^d{d,p) > such that 

lPMp(|Woo(0)|>n)<e-^^". 

Proof. As in the proof of Theorem 4.2, we stabilize rj in nested boxes 
Bk, and write [see (16)] 

IPm(|Woo(0)| >n)= lim P^(|Wfe(0)| > n) 

fc— >oo 

(oo 
J2 IP^drnWI = m) +P^(|Wfc(0)| = oo) 
— 
m=n 

The cluster Wfc(O) consists of the following types of sites: sites that have 
toppled, sites that did not topple but received at least one grain, and sites 
that did not topple nor received grains but which were nonempty in rj. The 
first two types of sites we can only find in the box -B^+ii but the third type 
we can also find outside this box. Outside the box B^+i, the configuration 
did not change yet, so restricted to Z'^ \ B^+i, we just have independent 
site percolation of nonempty sites. We take p so small that the density of 
nonempty sites is below the critical value for independent site percolation, 
so that for every k, |>Vfc(0)| is finite a.s. We write 

oo 

Pm(|Woo(0)| > n) = hm ^ F^(|W,.(0)| = m) 

m=n 

oo 

= 1™ E E iP/.(Wfc(o)=c), 

fe— +00 — , ~ — 

m=n\c\=m 

oeC 

and again derive a bound for P^(>Vfc(0) = C) using that, on the event {>Vfc(0) = 
C}, there must have been a certain minimal number of sand grains in C be- 
fore stabilization. Suppose Wfc(O) = C. If C contains a cluster of size mj > 1 of 
toppled sites, with > 2d boundary sites, then the number of grains in this 
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region of sites — after as well as before toppling — is at least rrif — 1 + mf, > 2d, 
so that the density in this region is at least ^M^- ^ might contain several 
of these regions, as well as nonempty sites that did not topple nor receive 
any grains. Thus, we cannot conclude more than that the density in C be- 
fore toppling was at least , which for p < corresponds to a large 

deviation of :^J2x&cV{x)- 

The rest of the proof proceeds the same as for Theorem 4.2. Note that 
the fact that we now sum m from n to co instead of to {2k + 1)'^ makes no 
difference for the outcome. □ 

Remark 4.6. For d=l, the critical density of nonempty sites is 1, but 
for all p we have that P^(?7(0) = !)<!. Therefore, Theorem 4.5 is valid for 
/3 < |. However, in d = 1 it is not hard to see that, for all p < 1, |Woo(0)| is 
finite a.s. Indeed, it is not hard to see that there ct.S. IS 3i positive density 
of pairs of neighboring sites which never topple (compare the last part of 
Theorem 2.8). This implies that |Woo(0)| is finite. 

5. Some open problems. It is clear that our results are a first step in 
the study of sandpile percolation and that many challenging open problems 
remain. We mention some of them: 

1. Infinite sandpile percolation clusters. Starting from a product measure 
Pp with density p on X, we know that, for p > 2d — 1, the measure is 
not stabilizable, and all sites will topple infinitely many times. For small 
p we have sandpile percolation clusters that look like subcritical clusters 
(of ordinary percolation). Are there values of p <2d — 1 such that there 
is an infinite cluster of toppled sites, but pp is still stabilizable? A guess 
would be that this happens as soon as p > d, the density of minimally 
recurrent configurations. If such a percolation transition occurs, is there 
a unique infinite cluster? 

2. Critical value. Let be a one-parameter family of Poisson product mea- 
sures with parameter p. Is pd-M.) <2d—l? Since it is conjectured that pc 
is the expected height in the critical sandpile model, in the limit of large 
volumes, this inequality should hold. However, we expect that, especially 
in high d, it will be easier to show a product measure with expected 
height "close to 2d — 1," and not concentrating on stable configurations, 
is not stabilizable. 

3. Behavior at the critical value. In d = 1 we proved that the Poisson mea- 
sure pp is not stabilizable at the critical value p=l. Is it true that, for 
all d > 1, pp is not stabilizable at the critical value p = pc, where pc is 
the critical value of the family of Poisson measures? 
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